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1 N-body simulations with 
17 million particles 

Dark Matter (DM) halos carry 
important information about the 
formation of galaxies and of many 
other extragalactic objects. In 
particular, the density profile and 
the clustering of DM halos are two 
very important ingredients for un- 
derstanding many observations of 
galaxies. In this talk, I present our 
new results for these two quanti- 
ties from a large set of high qual- 
ity N-body simulations. 

The simulations were gener- 
ated with a vectorized P 3 M (i.e. 
Particle Particle Particle Mesh) 
code on the the supercomputer 
VPP300/16R at the National As- 
tronomical Observatory of Japan. 
Each simulation uses 256 3 simu- 
lation particles. The code adopts 
the standard P 3 M algorithm. 
Twenty one simulations are gener- 
ated for three representative cold 
dark matter (CDM) models of 
galaxy formation, and eighteen 
for six scale-free cosmologies. The 
CDM models are specified com- 
pletely as regards the DM distri- 
bution by the density parameter 
Qo, the cosmological constant Ao, 
the shape P and the normalization 
as of the linear power spectrum. 
The scale-free models assume an 
Einstein-de Sitter universe (i.e. 
Qo = 1 and Ao = 0) and a power- 
law P(k) oc ft™ for the linear den- 
sity power spectrum. Table 1 sum- 
marizes the physical and simu- 
lation parameters used for these 
simulations. 

This set of simulations has cov- 
ered a large parameter space, 
which is often important in many 
studies. Each of our individ- 



ual simulations has similar mass 
and force resolutions to that the 
Virgo Consortium (Gus August 
in these proceedings) has re- 
cently obtained (except for their 
Hubble-volume simulation which 
has many more particles but still 
poorer absolute mass resolution 
than the CDM simulations pre- 
sented here), but our sample is 
much larger. The CDM simula- 
tions have been used to study the 
implications of the strong cluster- 
ing discovered by Steidel et al. 
(1998) of the Lyman Break galax- 
ies at high redshift (Jing & Suto 
1998). Many works based on these 
high quality simulations are be- 
ing done. As an incomplete list, 
we are 1) studying the DM dis- 
tribution in the strongly cluster- 
ing regime; 2) testing the Press- 
Schechter extension theories; 3) 
studying redshift distortion effect 
and investigating its potential to 
measure the cosmological param- 
eters; 4) analyzing internal struc- 
tures of DM halos; 5) investigat- 
ing galaxy formation by hand- 
inputing gas physics; 6) applying 
the simulations to real observa- 
tions. 



2 The density profile of the 
dark matter halos 

Recently Navarro, Frenk, & 
White (1996,1998; NFW) have 
found that the density profile p(r) 
of the DM halos follows a univer- 
sal form: 
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Table 1. List of simulations 

1. All simulations use 256 3 particles, except the low density models with 
power-law P(k) LSF1 to LSF4 use 200 3 particles; 

2. Particle mass m v in CDM models is units of h~ x Mq; 

3. Box size L in CDM models is units of h~~ 1 Mpc; 

4- The rightmost column is the number of realizations (samples); 



Model 


Ho 


Ao 
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0"8 


m p 


L 


Num. 


CDM1 


1.0 


0.0 


0.5 


0.6 


1.7 x 10 10 


100 


3 


CDM2 


0.3 


0.0 


0.25 


1.0 


5.0 x 10 9 


100 


3 


CDM3 


0.3 


0.7 


0.20 


1.0 


5.0 x 10 9 


100 


3 


CDM4 


1.0 


0.0 


0.5 


0.6 


4.6 x 10 11 


300 


4 


CDM5 


0.3 


0.0 


0.25 


1.0 


1.3 x 10 11 


300 
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CDM6 


0.3 


0.7 


0.20 


1.0 


1.3 x 10 11 


300 


4 


Model 


n 


Ao 


n 








Num. 


SF1 


1.0 


0.0 


1.0 
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SF2 


1.0 


0.0 


0.0 
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SF3 


1.0 


0.0 


-0.5 
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SF4 


1.0 


0.0 


-1.0 
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SF5 


1.0 


0.0 


-1.5 








3 


SF6 


1.0 


0.0 


-2.0 
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LSF1 


0.1 


0.9 


-1.0 
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LSF2 


0.1 


0.0 


-1.0 








2 


LSF3 


0.1 


0.9 


-2.0 








2 


LSF4 


0.1 


0.0 


-2.0 
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where r is in units of the ra- 
dius 7?200 within which the over- 
density around the halo center 
is 200, and the density profile 
depends only on the concentra- 
tion parameter c. The parame- 
ter c can be predicted for each 
cosmogonic model with a recipe 
given by NFW, and the predicted 
cnfw is a function of the halo 
mass only for a given cosmogo- 
nic model (i.e. the power spec- 
trum and the background cos- 
mology). These results are indeed 
very important and have many 
interesting applications. However, 
it is important to point out that 
NFW studied the density profile 



only for the halos which look in 
equilibrium, so it is unclear how 
much fraction of the DM halos 
really follows the universal form, 
since substructure is a common 
phenomenon for DM halos in all 
viable cosmological models (e.g. 
Jing et al. 1995). Even for ha- 
los which can be reasonably fit by 
the NFW profile, a dispersion in c 
seems inevitable because of a dif- 
ferent formation history. Quanti- 
fying this dispersion or more com- 
pletely the Probability Distribu- 
tion Function (PDF) of c is of con- 
siderable importance when apply- 
ing the NFW profile to interpret 
various observational facts. To an- 
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Figure 1. The probability distribution function of the concentration pa- 
rameter c. The right solid histogram is for the halos with dvimax < 0.15, 
the middle dotted one for 0.15 < dvimax < 0.35, and the left solid one for 
dvimax > 0.35. For each subset of the halos, the distribution of c can be 
nicely fitted by a lognormal distribution. 



swer these important questions, 
many halos with good resolution 
(certainly much more halos than 
that NFW had) are needed. Our 
simulations suit very well for this 
purpose since typically one model 
(either CDM or SF) has 300 to 
500 massive halos with at least 
10,000 particles within the virial 
radius. Both the force and the 
mass resolutions of our simula- 
tions are comparable to the NFW 
halos. 

We have selected DM halos us- 
ing the Friends-of-Friends method 
with the linking length 0.2 times 
of the mean particle separation. 
Then we choose the local poten- 
tial minimum as the center of 
the halo, and compute the den- 
sity p(r) in shells with logarith- 
mic thickness log 10 Ar — 0.1 from 
i?200 inward to rj, the force resolu- 
tion limit. The density profiles are 
then fitted with Eq. (1) to get the 
parameter c. Although we have 
computed the density profiles for 
the massive halos in all simula- 
tions, for the convenience of the 
discussion in this very limited 
space (and talk time), I present 
here the results only for one spe- 
cific model, the scale-free model 



of n = — 1. A detailed discussion 
of all the models will be given in 
a forthcoming paper (Jing & Suto 
1998, in preparation). 

Firstly we note that the good- 
ness of the fit varies strongly 
from halo to halo. Since the 
Poisson noise in the simula- 
tion density profile is negligi- 
ble, we define the goodness as 
the maximum relative deviation 
of the simulation p(r) from the 
fit pNFwir) in all the radial 
bins, i.e. dvimax = max{|(p(r;) — 
P NFW {n))/p NFW {r l )\}. About 
35, 50 and 15 percent of the ha- 
los have the maximum deviations 

dvimax > 35%, 15% < dvimax < 

35%, and dvimax < 15% respec- 
tively. With this classification, it 
is fair to say that about 35% halos 
could not be fitted by the NFW 
profile, because of too significant 
substructures; another 50% halos 
with less substructures can be rea- 
sonably described by the NFW 
profile, and the rest about 15% 
halos with the least substructures 
can be fitted by the NFW profile 
very nicely. 

The PDF of the fitted c is pre- 
sented in Figure 1. We choose 
c/cnfw as the abscissa in order 
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to correct for the mass depen- 
dence of the parameter c, though 
this correction for the mass range 
covered by the halos here is actu- 
ally tiny. The PDF is shown sepa- 
rately for the halos with different 
amount of substructures. For each 
subset of the halos, the PDF can 
be fitted by a lognormal function 



p(c)dc = 



1 



x exp - 



27TCT 

(lnc-lnc) 2 
2a 2 



dine. 



(2) 



As we discussed in the last para- 
graph, because of the poor fit 
for the halos with significant sub- 
structures, the PDF given by 
Fig. 1 for this subset is proba- 
bly not much meaningful and we 
shall not discuss it. But for the 
other halos which can be reason- 
ably fitted by the NFW profile, 
the PDF is indeed a very impor- 
tant quantity. Both for the most 
virliazed halos {dvi max < 15%) 
and for the halos with 15% < 
dvi m ax < 35%, the dispersion a 
of the PDF is 0.27. The mean 
value c is (0.93 ± 0.03)cnfw for 
the halos with dvi max < 15% and 
(0.84 ± Qm)c NFW for the halos 
with 15% < dvi max < 35%. The 
PDF of the parameter c depends 
weakly on the cosmological pa- 
rameters, and the above results 
can be applied to other cosmolog- 
ical models with flo >0.25. 

Our results are generally con- 
sistent with the results of NFW. 
Variation of the density profiles 
was also found by Kravtsov et 
al. (1998) but was discussed in 
a completely different way. What 
is new from our study is that 
we have quantified how much ha- 
los can really be fitted by the 
NFW profile and that we have de- 
rived the important quantity, the 
Probability Distribution Function 
of the concentration parameter c, 
for the halos which can be fit- 
ted by the NFW profile. These 
results would be important for 
properly interpreting many cos- 



mological observations which are 
closely related to the density pro- 
file of halos. 

Before concluding this section, 
it would be important to point 
out that because of the still lim- 
ited resolution of out simulations, 
we could not address the impor- 
tant problem that density pro- 
file might be much steeper than 
r _1 in the very inner region of 
the halos ( Fukushige & Makino 
1997; Moore et al. 1998). To def- 
initely answer this problem, one 
needs to simulate many halos with 
much higher resolutions of both 
the mass and the force than the 
present work. Our work in this di- 
rection is also in progress. 



3 The bias parameter of the 
dark matter halos 

The clustering of the DM ha- 
los is the basic block for model- 
ing and understanding the clus- 
tering of galaxies and of other 
extragalactic objects. We have 
recently measured the two-point 
correlation function of the DM 
halos for the scale free simulations 
with n = —0.5, —1.0, —1.5, and 
-2.0 (Jing 1998). The results con- 
firmed with unprecedented accu- 
racy that the bias of the halos b 
is linear in the linearly cluster- 
ing regime. For a scale- free model, 
the bias parameter b is then ex- 
pected to depend on the scaled 
mass M/M-,, only, where M* is 
a characteristic non-linear mass. 
In Fig. 2, I present the 6(M/M*) 
for different simulation outputs of 
each model. The excellent scal- 
ing exhibited by the simulation 
data assures that any numeri- 
cal artifacts indeed have negligi- 
ble effect on the results of Fig. 2. 
The simulation results agree with 
the analytical formula of Mo & 
White (1996) only for massive ha- 
los with M/M* > 1, but are 
significantly higher for much less 
massive halos. We also noted that 
there might be systematic differ- 
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Figure 2. The b 2 (M/M t ) for different evolution outputs (different symbols). 
The dotted lines are the analytical prediction of Mo & White (1996), and the 
solid lines are our fitting formula Eq. 3 which can accurately fit the simulation 
results. 



ence between the analytical pre- 
diction and our simulation results 
even at M/M* > 1 (see Fig. 2), 
despite not with a high statistical 
confidence. Our simulation results 
could be fitted very well by 

b{M ) = ( ^| + 1) (0.06-0.02n) 

xd + ^); 

v = (M/M*) n+3/6 , (3) 

where S c — 1.68. This fitting for- 
mula could be equally well ap- 
plied to the CDM models with 
only a very simple modification 
(Jing 1998). This finding has pro- 
found implications for the clus- 
tering studies of galaxies. Fur- 
ther discussions are given in Jing 
(1998). 
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